Experimental observations of dynamic critical phenomena in a lipid membrane 
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Near a critical point, the time scale of thermally-induced fluctuations diverges in a manner deter- 
mined by the dynamic universality class. Experiments have verified predicted 3D dynamic critical 
exponents in many systems, but similar experiments in 2D have been lacking for the case of con- 
served order parameter. Here we analyze time-dependent correlation functions of a quasi-2D lipid 
bilayer in water to show that its critical dynamics agree with a recently predicted universality class. 
In particular, the effective dynamic exponent Zctt crosses over from ~ 2 to ~ 3 as the correlation 
length of fluctuations exceeds a hydrodynamic length set by the membrane and bulk viscosities. 
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Lipids self-assemble in water to form sheets that are 
two molecules thick, within which the lipids are free to 
diffuse. When composed of several lipid species these 
two-dimensional (2D) liquid membranes can demix into 
coexisting liquid phases, termed Lq and Ld, over a range 
of temperatures and compositions, and can exhibit criti- 
cal behavior [iH^- Among 2D critical phenomena, com- 
position fluctuations in membranes are rather unique in 
that their large sizes and long decay times are accessible 
to optical microscopy. For example. Fig. 1 and supple- 
mentary movies show a vesicle (a spherical membrane 
shell) in which correlated regions reaching 10 /im persist 
for seconds [5]. Direct visualization of these equilibrium 
fluctuations has recently been used to show that static 
critical exponents for lipid membranes are consistent with 
the 2D Ising universality class [2 E] . Here we exploit the 
ability to visualize dynamics of these fluctuations to ex- 
amine for the first time the dynamic critical phenomena 
in this system. We find that although the statics are 2D 
phenomena, the critical dynamics are modified by hydro- 
dynamic coupling to the surrounding 3D fluid. 

Static critical exponents, which describe how observ- 
ables such as correlation length vary as the critical point 
is approached, are identical for all systems in a given uni- 
versality class, independent of their detailed microscopic 
physics [51 [7]. For example, although membranes have 
a conserved order parameter and ferromagnets do not, 
membranes exhibit static exponents v — 1.2 ± 0.2 and 
/3 = 0.124 ± 0.03, consistent with the expected 2D Ising 
values of = 1 and f3 = 1/8 Results in plasma 
membrane vesicles are also consistent with 2D Ising ex- 
ponents V = 1 and 7 = 7/4 [B]. Systems that are in 
the same static universality class can fall into different 
dynamic universality sub-classes determined by conser- 
vation laws constraining how fluctuations dissipate [S]. 
The critical exponent z for each dynamic subclass quan- 
titatively describes the scaling of the dynamics. It re- 
lates how the correlation time diverges as tempera- 
ture T approaches the critical temperature Tc, such that 
Ts (X \{T — Tc)/Tc\~'^^' where ly is the static critical ex- 
ponent. Experiments measure an effective exponent z^s 



that approaches z as T ^ Tc and ^ — oo. Dynamic 
sub-classes relevant to 2D systems with conserved order 
parameter are notable equally for their wealth of theo- 
retical predictions [5HTT| and for the lack of experiments 
that systematically test those predictions. 

Only a few previous measurements of dynamic critical 
exponents in 2D systems exist. Most experiments have 
been conducted on magnetic fllms. Using ferromagnetic 
fllms of ^ two monolayers, Dunlavy and Venus found 
i/z = 2.09 ± 0.06, with i/ ~ 1 [12]. Fewer experiments 
have been conducted on systems with conserved order 
parameter. Careful attempts to measure z were made in 
thin fllms of lutidine and water, but were unable to reach 
the 2D critical regime [13' . In plasma membrane vesicles 
from living rat basophil leukemia cells, fluctuation decay 
times were reported to be consistent with z « 2 j6j. 

Here we obtain z^s as T approaches Tc in a lipid mem- 
brane surrounded by water and compare to theory re- 
cently developed for an analogous system: a 2D critical 
binary fluid embedded in a non-critical bulk fluid [TUllllj . 




FIG. 1: Fluorescence micrographs of vesicles of diameter 200 
/im. (A) As temperature changes from T > Tc {T — 31.25°C, 
Tc ^ 30.9) to r ~ Tc (T = 31.0°C) fluctuations in lipid 
composition grow. Below Tc, at T = 28°C, domains appear. 
Scale bar — 10 ^m. (B) A movie of composition fluctua- 
tions within a vesicle above Tc. Large fluctuations persist for 
seconds (white arrows), whereas small ones disappear by the 
next frame (black arrow). Scale bar = 20 /im. 
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This new theory incorporates three essential features of 
Hpid bilayer dynamics: conserved order parameter, col- 
lective hydrodynamics, and hydrodynamic coupling be- 
tween the bilayer and bulk [101 HI]- Inclusion of only 
the first feature within an Ising model yields Model B, 
in which composition fluctuations dissipate through dif- 
fusion of microscopic constituents [H]. 2D Model B pre- 
dicts z = 4— 2/3 = 3. 75 [9], and numerical schemes give 
z = 3.80 and z ~ 3.95 [T3J [TS]. Inclusion of the first 
two features, such that collective hydrodynamic motion 
replaces single particle diffusion as the dominant mecha- 
nism of order parameter relaxation, yields Model H. 2D 
Model H with coupling to only 2D momentum modes 
predicts z « 2 [9]. Inclusion of all three features yields 
Model HC, where HC denotes hydrodynamic coupling of 
the membrane to the bulk. This new version extends 
Model H to account for modes in both the 2D membrane 
and the 3D bulk fluid, with the result that ^ =^ 3 [TOlITT] . 

Intuition for the role of the coupling between the mem- 
brane and bulk within Model HC can be gleaned from an 
approximation for 3D Model H by Kawasaki [TB] . Critical 
fluctuations are treated as spherical inclusions of diame- 
ter ^ that diffuse a distance ^ to equilibrate [9l H m [T6l [T7] . 
As such, correlation time varies as r ^ where 
£'(^) is the inclusion's diffusion constant in a non-critical 
fluid. In 3D, D(t) ^ where r is the inclusion's ra- 
dius. Using T cx \[T -T^)IT^\~'''' cx yields z w 3. 
A more sophisticated theoretical treatment gives z ~ 
3.065 p^. Applying the same reasoning to 2D Model 
H, in which diffusion of inclusions has only a logarithmic 
dependence on r, yields z « 2. Again, more sophisti- 
cated treatments produce similar values; see [S] for more 
detail. This argument can be extended to predict the 
value z should take in a 2D critical system embedded in 
a bulk fluid. Classic work by Saffman and Delbriick ex- 
amined diffusion of an inclusion in a 2D liquid of viscosity 
?72D immersed in a bulk fluid of 3D viscosity rysD, where 
hydrodynamic length Lh = '72D / ?y3D is an important pa- 
rameter [19l[20]. When ?■ ^ Lh, dissipation is primarily 
into the bulk and D(r) c>c I/r as in 3D Model H. When 
r <C ih, dissipation is primarily into 2D hydrodynamic 
modes and D(r) oc ln(Lh/r'), similar to 2D Model H. 
Two groups have independently noted that when Lh is 
considered, for a 2D critical binary fluid embedded 
in bulk liquid crosses over from ZoS ~ 2 when ^ ^ to 
Zeff ~ 3 when ^ > Lh CHI HI] ■ 

The next four paragraphs demonstrate that experi- 
mental results here are in excellent agreement with the 
recent predictions of Model HC, namely that Zg// crosses 
over from '--^ 2 to ~ 3 as T — )■ Tc and ^ — )• oo. Further 
experimental details follow the results. 

A time series of the order parameter, 'm{r,t), was 
extracted from videos of vesicles collected via fluores- 
cence microscopy. For membranes, m{r) is the deviation 
from average composition as reported by an image's pixel 
grey scales. A time-correlation function C(r, r), and its 



Fourier transform in space, the structure factor S{k,T), 
were calculated for each wavenumber k. 

Curves of S{k,T)/ S{k,0) vs k^^'W were plotted for 
a range of Zcs values. Fig. 2B illustrates how the 
correct z^s was identified: for a single value of Zcg, 
all experimentally-measured curves at different k values 
(Fig. 2A) collapsed most fully onto a single curve, here at 
Zoff = 2.8±0.2. Fig. 3A shows Zes values extracted in this 
manner from data over the entire measurable range of 
correlation lengths. In Fig. 3A, Zeg rises from from near 
2 to near 3 as T ^- Tc, in accord with Model HC [TU l fTT ] . 

Fig. 2C-D validates this method by showing that stan- 
dard simulations of Model B Kawaski dynamics that are 
blurred to mimic experimental limitations and then an- 
alyzed in the same way as the experimental data give 
z — 3.6 ± 0.2 in agreement with the expected value of 
z = 3.75 (see |S] for details). Simulations were run on 
a 400x400 bi-periodic square lattice. Blur was achieved 
by averaging snapshots over 200 consecutive Monte Carlo 
sweeps, leaving a break of 800 sweeps without snapshots, 
and repeating the process, which reproduced the effects 
of a camera shutter opening for 100 ms of every 500 ms. 

Excellent agreement between predicted and measured 
structure factors provides even stronger evidence that 
Model HC describes critical dynamics in membranes. In- 
aura and Fujitani give a prediction for the entire time- 
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FIG. 2: (color online). (A and B) Rescaling experimental data 
closest to Tc by k^r collapses all curves to Zea = 2.8, consistent 
with Model HC. Normalized structure factors are shown for 
^ = 13 ± 2.2 /im and three video rates: 10 frames per second 
(fps, solid lines), 2 fps (short dash), and 0.5 fps (long dash). 
Colors denote wavenumbers k = 1.1 fim~^ (top curve, blue) 
to 3.0 fim~^ (bottom, red). (C and D) Simulations solely 
to verify technique. Structure factors of Kawasaki dynamics 
at T — Tc blurred in time to mimic experimental limitations 
collapse at Zea ~ 3.6 ± 0.2, consistent with z — 3.75 for 2D 
Model B. Colors range from k — 1.1 fj,m~^ to 3.1 /im"^. Insets 
show collapses used to determine bounds for ZcB and failure 
of collapse at Zeff = 3. 
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dependent structure factor S{k,T) for Model HC, tak- 
ing as input 772D, and a mean- field approxima- 
tion for the static structure factor, S{k,0) [11]. The 
ratio S{k,T)/ S{k,0) and its decay time will be com- 
pared between theory and experiment below. A fea- 
ture of S{k,T) / S{k,0) is that it needs no correction 
due to the microscope's point spread function. Ratios 
of S{k,T)/S{k,0) in the critical Ising model and in the 
mean-field approximation are similar, and do not depend 
strongly on correlation length, as will be shown in a fu- 
ture manuscript. 

The HC model with — 6 fim fits the data over all 
experimentally accessible wavenumbers. Fig. 3C shows 
the ratio S{k,T)/ S{k,0) at wavenumbers 1.1 /xm~^ and 
3.8 ^m~^. Fig. 3B shows decay times, defined as when 
S{k,T)/S{k,0) = e~^. All other models are excluded. 
Fig. 3A rules out 2D Model H because the measured ZeS 
rises close to 3, well above the predicted value of 2 for 
2D model H. Fig. 3B-C rules out 2D Model B because 
measured decay times are orders of magnitude shorter 
than the model predicts. 

Developing Model HC to completely describe mem- 
branes requires determining only membrane viscosity. 
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FIG. 3: Data is in excellent agreement with Model HC. 
(A) Filled symbols: Dynamic exponent ZcB from scaling col- 
lapse of experimental data as in Fig. 2A-B. Open symbols: 
Model B simulation in which ZeS approaches ~ 3.75. (B) 
Decay time, defined as when S{k,T) / S{k,0) = e~^. Large 
symbols indicate wavenumbers 1.1 and 3.3 fim~^ . (C) Nor- 
malized structure factors r)/S(fc, 0). In panels B and C, 
experimental data is denoted by symbols, 2D Model B by a 
grey line. Model HC (HC) with Lh = 6 /xm by a solid line and 
Model HC with Lh = by a dashed line. 



rj2u, as an input since the viscosity of water, 7730, is 
known. Fig. 3C indicates that 772D must be nonzero. In 
the small k limit, Inaura and Fujitani [11] predict a struc- 
ture factor that depends only on 7730- This parameter- 
free prediction, equivalent to taking 772D = 0, underesti- 
mates time decays by a factor of 5 — 10 (dashed curve. 
Fig. 3C). Setting the unknown 7720 (or equivalently, Lh 
— ??2d/??3d) as a single fit parameter within Model HC 
over the entire measured range of k yields L^ = 6.0 ± 
1.5 fim. This value is within the range found by tracking 
diffusion of liquid domains across vesicle surfaces [HJ |22] 
and is similar to values (2-4 /xm) found by other methods, 
albeit for different lipid mixtures [23l [24] . An essentially 
equivalent method of finding Lh is to calculate the Model 
HC structure factor using the formalism of Hohenberg 
and Halperin [S], and to thereby extend Model HC to 
incorporate Ising rather than mean field statics. Within 
experimental uncertainty, this modest change has no ef- 
fect (Lh = 5.5 ± 1.5 /im). This and other extensions of 
Model HC, each leading to small corrections to the ratio 
5(fc, r)/5(fc, 0), will appear in a future manuscript. 

Using lipid bilayers to measure critical exponents in- 
troduces both complexities and advantages, which are 
outlined further in [S]. The first complexity is that the 
simplest bilayers that exhibit critical phenomena contain 
ternary lipid compositions. Strictly speaking, the ternary 
mixtures used here pass through isothermal critical mix- 
ing (plait) points rather than critical (col) points. A fea- 
ture of 2D systems is that, unlike in 3D systems, no mea- 
surable change in critical exponents arises from the pres- 
ence of a third component. Briefly, a small correction to 
scaling arises in systems that contain a third component 
at flxed composition rather than fixed chemical potential. 
Hence, Tc changes, and many effective critical exponents 
are renormalized by a factor of 1/(1 — a), as discovered 
by Widom [2] and generalized by Fisher [3]. For the 2D 
Ising case here, where a = 0, theory predicts only a log- 
arithmic correction to singular behavior [2]. The second 
complexity is that when T is changed (as required in pre- 
vious studies to find v and /3 [3| , but not required here) , a 
bilayer with fixed composition does not necessarily follow 
a path with constant (m(r)) |27j . However, since mem- 
brane phase diagrams are relatively symmetric over the 
range of temperatures probed and since measured values 
of I' and /3 were consistent with the 2D Ising model [5], 
deviations from a path of constant {m(r)) are likely mi- 
nor. 

The first advantage of using lipid bilayers is that it 
avoids challenges of other systems. For example, lipid 
monolayers have confounding effects of dipole interac- 
tions, and the task of achieving simultaneous tunability 
and stability of surface pressure in a stationary mono- 
layer is formidable. The second is that correlation 
lengths are large, partly because ^0 within the relation 
^ = ^0 |(T - Tc)/Tc\~'' is on the order of the length of a 
lipid molecule rather than of an atom. Separately, there 



4 



is an advantage in using 2D (or quasi-2D) experimental 
systems over 3D systems. The critical region is larger 
in 2D liquid-liquid critical systems than in analogous 3D 
ones, partially due to differences between critical expo- 
nents in 2D vs 3D Ising classes {1/ = 1 and /3 = 1/8 in 
2D vs sa 0.630 and l3 sa 0.325 in 3D [5]). 

Methods used to produce the results in Fig. 3 fol- 
low. To optimize movie quality, vesicles were spherical, 
free-floating, unilamellar, of radius >100 ^m, and electro- 
formed by standard methods detailed in j3] . Vesicles were 
formed from mixtures along a line of plait points centered 
at 30% diphytanoylphosphatidylcholine (DiPhyPC), 20% 
dipalmitoylphosphatidylcholine (DPPC) and 50% choles- 
terol (chol), with 0.5% fluorescent dye Texas red di- 
palmitoylphosphatidylethanolamine (TR-DPPE) . Only 
vesicles near a plait point were analyzed, identified by 
micron-scale composition fluctuations visible over the 
largest observed range of temperatures (> 1°C) and by 
equal areas of coexisting liquid phases below Tc. Each 
vesicle analyzed fell on a slightly different plait point, so 
each had a slightly different Tc [28] . 

Images of membranes were captured via an epifluo- 
rescence microscope with a temperature-controlled stage 
and a mercury lamp source. Light exposure was mini- 
mized by employing a SmartShutter (Sutter Instrument, 
Novato CA) controlled through NIS-Elements (Nikon, 
Melville NY) and by recording movies for at least two 
different frame rates at each temperature. Each frame 
was exposed 100-150 ms, with the shutter open 10 ms 
before and after exposures. Movies were collected from 
high to low temperature in steps of ^ 0.2°C, equilibrated 
for at least 2 min. No consistent trend in intensity was 
observed throughout each movie, implying that the low 
light procedures used here eliminated signiflcant photo- 
bleaching. To correct for lamp flickering, mean brightness 
was subtracted from each frame. Spatial intensity gra- 
dients due to other vesicles outside the focal plane were 
removed by a long wavelength filter of 100 pixels. 

Images were analyzed via custom MATLAB code (The 
Mathworks, Natick, MA) . Vesicles were tracked and cen- 
tered to remove drift (typically < 25/im/min). By eye, 
features exhibit no net translation, which implies no sig- 
nificant vesicle rolling. No difference in mean intensity or 
noise between pixels at edges vs. centers of cropped im- 
ages was observed, implying that vesicles are so large that 
membrane curvature over images can be neglected [3]. 
Curvature corrections in smaller vesicles were minor [S]. 

The structure factor S{k,T), the Fourier transform in 
space of the time-dependent correlation function, was 
found as previously described [3 [H]- Briefly, a dis- 
crete transform was performed for each movie image, 
with a buffer of zero values to correct for image non- 
periodicity. Transformed images were divided by the mi- 
croscope's finite point spread function to yield m{k,t). 
The dynamic structure factor was generated at each r 



by S{k,T) = 1/2 (rn{k,t)m{k,t ± t)^ , where m{k,t) is 

the complex conjugate of m{k,t) [30] ■ S{k,T) was then 
radially averaged to yield S{k,T). 

Structure factors were employed in two ways. First, 
correlation lengths, ^, were found by analyzing structure 
factors at r = 0. Specifically, a one-parameter fit for ^ 
was made until all data for k^'^^^^S{k) vs. fc^ collapsed 
onto the single curve for the exact numerical solution of 
the 2D Ising model [31 [3T]. Second, effective dynamic 
scaling exponents, ZeS, were found by collapsing curves 
of S{k, t) (see results above and [S] for details). Collapse 
works because, according to the dynamic scaling hypoth- 
esis, structure factors within the scaling regime can be 
written in the form 5(fc,r,C) = fc-2+'?f7((fc^)-i^ /fev) 
where O is a universal function of (fc^)"^ and fe^r [3]. 
Near Tc, where {k£^)~^ is near 0, curves of S{k, T)/S{k, 0) 
vs k^T collected over many wavenumbers k should col- 
lapse via a one-parameter fit to produce the correct 
value of z. Here, il can also depend on /cLh, so that 
5(fc,T,0 = k-^+in{{k^)-^ , k'^T, kLi,). For collapses in 
Fig. 2A-B, Zoff refers to an effective z value which varies 
as C/ih is changed. In Fig. 3B-C, comparing the entire 
form of the structure factor to theoretical predictions di- 
rectly verifies the value of z as well as the dependence of 
the universal function on /cLh and k^r. 

Summary: Directly imaging composition fiuctuations 
enables measurement of effective dynamic critical ex- 
ponents of a lipid membrane embedded in bulk wa- 
ter. Experimental structure factors are in excellent 
agreement with an emerging theoretical prediction in 
which 3D hydrodynamics affects critical slowing down 
in a 2D membrane. The theory invokes hydrodynamic 
coupling between the membrane and bulk fluid such 
that Ising degrees of freedom are coupled to momentum 
modes [TUHH]. As predicted, a shift in z^s from ~ 2 to 
~ 3 as T — )■ Tc and ^ — 00 is observed. 
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SUPPLEMENTARY INFORMATION 
Introduction 

This supplemental section contains details regarding 
analysis, calculation, and theory not found in the main 
text. Section[]contains additional information on analysis 
procedures. Section[]explains why Fisher renormahzation 
does not affect the observed static critical exponents in 
2D Ising systems, and why the critical region is larger 
in 2D liquid-liquid membrane systems than it would be 
in analogous 3D ones. Section [] details the Ising model 
simulations that we used to verify our analysis. Section [] 
reviews previous theoretical results for purely 2D critical 
hydrodynamics . 

Details of Analysis 



In the definition of S{k, r) = yri{k, t)m{k, t + t)^ , the 

value of S{k, r = 0) is guaranteed to be real as an expec- 
tation value and also from run to run since each term 
appears with its complex conjugate. For t 7^ 0, the ex- 
pectation value of the imaginary part of S is guaranteed 
to be real by time-reversal invariance expected of sys- 
tems in equilibrium. However, each term contributes an 
imaginary component. Adding the complex conjugate 
and dividing by two leads to an effective measurement 
of l/2(5(fc, r) -t- S{k, — r)), which for an equilibrium sys- 
tem is equal to S{k,T). The inclusion of zeros required 
to pad the raw data in m{k,t) introduces a small error 
in the calculation of S{k,T), which can be corrected for 
by dividing S{k, r) by the correlation function of pure 
ones and zeros in the real space correlation function. For 
calculations in fc-space, there is no simple correction (the 
real space correction has an ill-posed Fourier transform 
and so introduces unacceptable noise in fc-space). Never- 
theless, any correction is expected to be small (values of 
S{k, t) were similar when calculated with vs without zero 
padding) . More importantly, any correction would cancel 
out of the main results presented here, where 5(fc, r) is 
divided by S{k,T = 0). For movies at the slowest frame 
rate, 0.5 fps, noise in iS'(fc,0) caused an offset from the 
rest of the structure factor. For calculations made with 
those data sets, the measured value of ^(fc, 0) was re- 
placed by the value extrapolated from an exponential fit 
to the 2nd through 5th points in S{k,T). 

Experimental advantages of lipid bilayers 

This section explains why Fisher renormahzation does 
not affect the observed static critical exponents in 2D 
Ising systems, and why the critical region is larger in 
2D liquid-liquid membrane systems than it would be in 
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analogous 3D ones. 

Widom-Fisher rescaling leads to only an immeasurably 
small correction to 2D critical exponents: Concentration 
fluctuations consistent with static 2D Ising critical expo- 
nents were previously observed in membranes over our 
entire range oi pj. As we show below, this observa- 
tion is not at odds with the fact that our ternary system 
is subject to rescalings first discovered by Widom |2] and 
generalized by Fisher [3 . Rescaling corrections apply 
to any system with a quantity whose chemical potential 
smoothly affects the critical temperature. When a com- 
ponent is instead held at fixed composition (as our three 
components are) then the observed critical behavior re- 
ceives non-analytic corrections, essentially because the 
chemical potential of the third component has singular 
behavior near the critical point when held at fixed com- 
position. As a result, the singular form of the coexistence 
curve near the critical point is changed from its usual ex- 
ponent P to /?' = — a). Here a is the static critical 
exponent for specific heat. Other critical exponents that 
relate singular behavior of a quantity to the distance in 
temperature from the fixed point (for example a, /?, v 
and 7) receive similar corrections. For example, the spe- 
cific heat exponent itself becomes a' = a/(l — a)[3]. The 
rescaling correction is not confined to ternary systems: 
binary systems held at fixed density, rather than fixed 
pressure undergo similar rescalings, as does any system 
in which a density variable is held fixed in a phase dia- 
gram rather than it's conjugate field. 

Fortunately for the current study, Widom-Fisher 
rescaling in the 2D Ising model leads to an immeasurably 
small change in the singular behavior and no change in 
the critical exponents themselves. As noted above, a is 
the static critical exponent for specific heat, C, which di- 
verges as C ~ ((T - Tc)/Tc)-°' . In the 2D Ising model, 
specific heat diverges as C ~ log((T — Tc)/Tc), which is 
slower than any power law divergence, so that a is said to 
be zero. As explicitly discussed in both [5] and there 
is potentially a logarithmic correction to the singular be- 
havior of quantities whose critical exponents are usually 
multiplied by (1 — a)~^. For example, the correlation 
length ^, which is usually written as ^ oc ( )~'^, be- 
comes ^ oc ( ~ iT-T^){T-T^) ^-i^Q^f^-g]^ rescaling. In this 
case, v' = V, so rescaling does not change the critical 
exponent. For a system described by the 3D Ising model 
(rather than the 2D Ising model as in the current study) , 
a K, 0.11 131 H], such that the effect of rescaling is small 
but observable in the critical exponents. 

Although the effect of rescaling on the dynamic ex- 
ponent z is not discussed explicitly in the literature, we 
expect that z would not be affected by rescaling even 
for systems in which a ^ 0. The dynamic exponent z 
describes scaling of the time scale as the length scale is 
changed, with cx The product vz describes crit- 
ical slowing with respect to temperature, where Tg cx 



{ '^TjT'" y . As such, vz does describe the singular be- 
havior of a quantity (here the time scale) as temperature 
is changed and it receives a correction of through 
the parameter v. 

The critical region is larger in 2D liquid-liquid mem- 
brane systems than in analogous 3D ones. The point 
here goes beyond considerations that as dimension de- 
creases, the critical region becomes larger due to changes 
in the Ginzburg temperature, which is the temperature 
at which a system crosses over from mean-field to crit- 
ical behavior [Ij. Correlation lengths of f ~ 10/zm are 
regularly observed in vesicle membranes. These corre- 
lation lengths are larger than those typically observed 
in 3D binary mixtures, even though control over lipid 
composition in membranes is coarser. This is because: 
(1) lipids are molecules with length scales of ^0 ~ Inm, 
whereas the atoms employed in many studies of 3D crit- 
ical phenomena are an order of magnitude smaller, with 

^ O.lnm, and (2) differences between critical expo- 
nents in the Ising classes are favorable to 2D systems. 
In 2D, 1 and /3 1/8, whereas in 3D, v w 0.630 
and P ~ 0.325 [5]). The scaling form for correlation 

length can be written as ^ = ^Qt~^U ^(|0| /0o)^^^ , 

where t — (T — Tc) /T^ is reduced temperature, is a 
molecular scale, and 'magnetization' is a function of 
composition normalized by 4>q, which is roughly the dif- 
ference in composition between the two low temperature 
phases far from the critical point. lA{x) is a universal 
function of its dimensionless argument x. lA{x) has a 
maximal value at a; = and decreases to as its argu- 
ment increases. Tuning a system's correlation length to 
^, requires tuning temperature to within approximately 
AT ~ Tc{£,/ ■ Consider a membrane with values 
observed here: ^ ^ 10 ^m, ^0 1 nm, and Tc ^ 300K. 
In 2D, AT ^ .03 K, which is experimentally achievable. 
For 3D systems in which is typically 1 A, a 10/xm 
correlation length would require tuning temperature to 
a much higher precision of AT ~ 10~^K. Similarly, ex- 
treme accuracy in composition is not required in 2D to 
observe fluctuations of ^ ^ 10 /im. For a membrane to 
lie in the critical region, composition must be tuned such 
that the argument oiU is of order 1. Given the difference 
between critical exponents in 2D vs 3D, far less precise 
experimental control of composition can be tolerated in 
2D than in 3D. 



Simulation details 

Simulation procedures were standard O [7] and briefly 
explained here. The standard Ising Hamiltonian given by 
H — — j} ^i^j '^^s used, with spin variables Si — ±1 
and summation over the four nearest neighbors (j) of ev- 
ery state (i). Temperatures were in terms of the exact 
critical temperature given by the Onsager solution |S], 
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Tc = 2/log(l + so that a reduced temperature 
t ~ {T — Tc)/Tc corresponds to a simulation tempera- 
ture of Ts„„ = 2.269(1 + i). In this section, T and AiJ 
correspond to temperature and to the change in energy 
between initial and final states, respectively. Both are in 
dimensionless units. In a Monte-Carlo 'sweep', 160,000 
(400^) pairs of spins were proposed to be swapped, such 
that each spin was proposed twice. Metropolis spin ex- 
changes were used; each pair was exchanged or not to 
satisfy detailed balance [H [7| . If the resulting configura- 
tion was lower in energy, the exchange was accepted. If 
energy increased, the exchange was accepted stochasti- 
cally with probability exp(— AiJ/T). 

Note that any dynamics that satisfy detailed balance 
will lead to the same equilibrium ensemble of configura- 
tions [7]. To rapidly equilibrate the system, 'nonlocal' 
moves were employed in which each of a pair of spins 
were chosen from all sites on the lattice. Equilibration is 
very rapid using these nonlocal dynamics since they ap- 
proximate "Model A" for large systems where z is near 
2 [5]. The system was equilibrated for 100, 000 sweeps us- 
ing nonlocal moves starting from a distribution that con- 
tained the desired fraction of up spins but was otherwise 
random. 100, 000 sweeps is much longer than the decay 
time of the slowest decaying system used here. The decay 
time is approximately 1000 sweeps at l.OSTcj which can 
be seen qualitatively by inspecting successive snapshots 
or quantitatively by inspecting the decay of time depen- 
dent correlation functions. Once the system was equili- 
brated, dynamics relevant for the locally conserved order 
parameter (Kawasaki Dynamics) were employed. In this 
case, a single spin and one of its four nearest neighbors 
were chosen to form a pair proposed to be swapped. 

Predictions for binary liquids in 2D 

Model H for binary fluids in 2D predicts z w 2 using 
z = 4 — 1] — xx, where 77 = 2/3 is a static critical exponent 



and x\ must be calculated from an epsilon expansion 
(where e = 4 — Z?, and D is the number of dimensions). 
This yields xx = 18/19(1 - (constant)e + Oe^) where the 
constant is either 0.033 [lUI or 0.039 [TT]. Since the con- 
stant is small, it is plausible that the expansion applies 
even when e = 2, yielding z — 2.00 (which also arises 
from a much simpler mean field argument) or 1.98. Sim- 
ulations in 2D binary liquids are reportedly challenging 
and we know of none that either verify or contradict the 
prediction that z = 2. Measurements in bulk 3D liquids 
far from find that z^s = 2 (e.g. [HHH])- 
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